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1. Introduction 

The spectral energy distribution (SED) of an ultraluminous infrared galaxy 
(ULIRG) at wavelengths between about 10 and 1000 /im is dominated by 
the thermal emission from interstellar dust, which is heated by absorbing 
blue and ultraviolet light from young stars and /or an active galactic nucleus 
(AGN). Mid- and far-IR radiation from galaxies at redshifts z <>1 has been 
detected directly by the IRAS and ISO space-borne telescopes. The red- 
shifted dust emission from more distant galaxies can be detected efficiently 
at longer submillimetre wavelengths. The only one with a known redshift 
and a constrained mid-IR SED is SMM J02399— 0136 at z = 2.8 (Ivison 
et al. 1998). High-z galaxies were also detected by IRAS, but only those 
with flux densities enhanced by gravitational lenses. Three such galaxies 
are known - IRAS F10214+4724 (Rowan-Robinson et al. 1991), H1413+117 
(Barvainis et al. 1995) and APM 08279+5255 (Irwin et al. 1998). The SEDs 
of the four high- z galaxies, and two well studied low-z ULIRGs - Arp 220 
and Markarian 231 (Klaas et al. 1997; Rigopoulou et al. 1996; Soifer et al. 
1999) - are shown in Fig. 1. 

2. Spectral energy distributions 

The difference in the slope of the mid-IR SEDs of Arp 220 and Mrk 231 is 
interpreted as evidence for an AGN in the core of Mrk 231. Powerful ionizing 
radiation near to the AGN would heat a very small fraction of the dust 
grains to high temperatures, boosting the flux density of the galaxy at short 
wavelengths and thus producing a shallower SED. There have been various 
attempts to investigate the mid-/far-IR SEDs of dust-obscured AGN using 
radiative transfer models (Granato et al. 1996; Green & Rowan-Robinson 
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Figure 1. The SEDs of two low-z (left panel) and four high-z (right panel) ULIRGs. Note 
that the mid-IR SEDs of APM 08279+5255 (Lewis et al. 1998; Downes et al. submitted), 
H1413+117 and IRAS F10214+4724 are probably modified by differential magnification. 
The SEDs are described by a power law, /„ oc v p , at short wavelengths and by a blackbody 
spectrum at temperature Td, modified by a dust emissivity e„ oc with /? ~ 1 at long 
wavelengths. The best fitting values of T d and p are listed in Table 1. 

1996), including the effects of an obscuring dust torus, the orientation of 
which determines the optical depth to the AGN, and thus its appearance 
in the optical (Antonucci 1993) and mid-infrared wavebands. 

However, given that a power-law spectrum convincingly accounts for 
the limited mid-IR data available, it is reasonable to avoid such details and 
assume that the SED of the galaxy can be described by the superposition 
of the emission from dust clouds at different temperatures, corresponding 
to different distances from the nucleus in the case that a significant fraction 
of dust heating in a galaxy is due to an AGN. If the dust temperature is 
assumed to depend on the distance from the nucleus r as T r (r) oc r v (for a 
blackbody in equilibrium with an unobscured point source r] = —0.5), and 
the mass of dust enclosed in the spherical shell between radii r and r + dr 
is m r (r) oc r 7 , then dust at each temperature can be associated with a 
(5-function spectrum, 5[v — vq\, in which vq ~ (3 + (5)kT/h (Blain & Longair 
1993), and the mid-IR SED, 

/ f max 
m v {r)T v {rf +p 5[v{l + z)-v Q ]dr. (1) 
min 

j3 is the index in the function that describes the spectral emissivity of 
dust e v oc (3 ~ 1-1.5. Dust temperatures between T m - m = 40 K and 
T ma ,x = 2000 K correspond to emission at wavelengths from about 90 to 
1.8 /urn. f u thus evaluates to a power law v v with p = 3 + (3 + [(7 + l)/r/]. 
More details of the calculation can be found in Blain (1999). 
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TABLE 1. The redshift, optical spectral type (Antonucci 1993), restframe dust tem- 
perature, and mid-IR spectral index p - where /„ oc v p - for the galaxies whose SEDs 
are plotted in Fig. 1. For comparison, the spectral index p ~ —2.0 to —2.5 for samples 
of low-redshift IRAS galaxies, either those selected at a wavelength of 25 /xm (Xu et 
al. 1998) or for the most luminous selected at 60 pm (Sanders & Mirabel 1996). 



Name 


Redshift z 


Type 


Restframe dust 


Mid-IR spectral 








temperature T d / K 


index p 


Arp 220 


0.02 


2 


50 ±1 


-3.6 ±0.2 


Markarian 231 


0.03 


1 


47 ± 1 


-1.9 ± 0.2 


IRASF10214+4724 


2.3 


2 


76 ±4 


-1.7 ±0.3 


H1413+117 


2.6 


1 


75 ±4 


-1.1 ±0.2 


SMM J02399-0136 


2.8 


2 


53 ±4 


-1.7 ±0.2 


APM 08279+5255 


3.9 


1 


107 ±8 


-1.0 ±0.3 



Two of the three high-z galaxies with IRAS detections - APM 08279 
and H1413 - have much flatter mid-IR SEDs as compared with a typi- 
cal low-z IRAS galaxy or SMMJ02399: see Fig. 1. Although SMM J02399 
lacks a detection by IRAS, its mid-IR SED is constrained by a 15-//m ISO 
measurement of L. Metcalfe (Ivison et al. 1998). Given the limited mid-IR 
data for SMMJ02399, it is certainly possible that its mid-IR SED could 
be steeper than the value of p = —1.7 listed in Table 1. Another high-z 
source, the brightest detected in a 850-/xm survey of the Hubble Deep Field 
(Hughes et al. 1998) has a reported 15-/um flux density upper limit that 
indicates a mid-IR SED with p < —1.7. Note that the SEDs of these galax- 
ies are unlikely to be affected by contaminating sources picked up in the 
different observing beams used at each wavelength. The beam sizes of the 
telescopes are about 5, 10, 25, 40, 7 and 14arcsec at 12, 25, 60, 100, 450 and 
850 fim respectively. Any contamination from other sources nearby on the 
sky would thus be most significant at wavelengths of 60 and 100 /xm, and 
lead to a artificial steepening of the mid-IR SED. The low-redshift sources 
listed in Table 1 are flagged as unresolved at these wavelengths in the IRAS 
catalogue, apart from Arp 220 at 12 //m, and so aperture photometry errors 
are not expected to modify their mid-IR SEDs significantly 

The shallow SEDs of APM 08279 and H1413 could be due to either an 
unusually large fraction of hot dust in these galaxies, or a relatively face-on 
aspect of an obscuring torus, as indicated by a type 1 optical spectrum, 
which would allow a direct view of the hot dust in the inner regions of 
the torus. However, both are known to be gravitationally lensed, and a 
systematically greater magnification for hotter dust components would also 
increase the flux density at shorter wavelengths and flatten the spectrum 
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(Eisenhardt et al. 1996; Lewis et al. 1998). This situation would arise very 
naturally if the hotter dust clouds are smaller and more central, as described 
above. 

For a large magnification to occur, a distant source must lie close to a 
caustic curve of a gravitational lens. The magnification is formally infinite 
on such a curve, but an upper limit A max is imposed to the magnification 
if the source has a finite size d (Peacock 1982); i max oc d^ 1 . This rela- 
tionship between size and magnification holds regardless of the geometry of 
the source. The diagnostic feature of such a situation is the production of 
multiple images of comparable brightness. This is more clearly the case for 
both APM 08279 and H1413 than for IRASF10214; however, IRASF10214 
lies very close to the tip of an astroid caustic (Broadhurst & Lehar 1995), 
and so differential lensing would still be expected to flatten its mid-IR SED. 

The effects of differential magnification on the mid-IR SED can be cal- 
culated by including a factor of A max (r) oc r _1 in the integrand of equation 
(1). In this case, the mid-IR spectral index p = 3 + (3 + (7 /rj), reduced by 
Ap = I/77 as compared with the unlensed value. Note that the form of the 
modification to the SED depends only on the value of 77, and thus the form 
of T{r). It is independent of the value of 7 and the form of m(r). Following 
the same approach, but assuming an exponential decrease of T(r), Ap ~ 1. 

3. Discussion 

The mid-IR SEDs of the high- and low-z AGNs SMM J02399 and Mrk231 
have spectral indices p ~ —1.8. Despite being lensed by a cluster of galaxies 
SMM J02399 does not lie close to a caustic, and so is not magnified differ- 
entially. The strongly lensed distant galaxies APM 08279 and H1413 have a 
flatter spectrum with p ~ — 1.1. If the difference is caused entirely by differ- 
ential magnification, rather than by a torus orientation effect (Granato et 
al. 1996), then, as Ap(= I/77) ~ —0.7, 77 ~ —1.4. Reasonably, this value of 
77 corresponds to a temperature profile that declines more rapidly as com- 
pared with that of an unscreened blackbody. To match the observed SEDs, 
7 ~ 7.3, indicating that most of the emitting dust is relatively cool. 

IRASF10214 has a steeper mid-IR SED, similar to those of Mrk231 
and SMMJ02399, and a much smaller optical flux density as compared 
with APM 08279 and H1413 (Lewis et al. 1998). These features can both 
be explained if the optical depth of dust extinction into the central regions 
of IRASF10214 is sufficiently large to obscure the AGN from view even 
in the mid-IR waveband, consistent with an inclined optically thick dust 
torus and the type-2 optical spectrum of this galaxy. Alternatively, the 
intrinsic unmagnified SED of IRASF10214 could be steeper than those of 
APM 08279 and H1413, just as these two galaxies could have intrinsically 
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flat mid-IR SEDs. 

Does a value of rj = — 1.4 correspond to a plausible size for the emitting 
region of a ULIRG? Dust sublimes at a temperature T ~ 2000 K, which 
is the equilibrium temperature of a blackbody heated by a 1O 13 -L0 point 
source at a distance of 0.6parsec. If rj = —1.4, then T falls to 100, 50 and 
30 K at distances of 5, 8 and 14parsec respectively. This is all in agreement 
with observations of low-z ULIRGs, in which the emitting region is less 
than several hundred parsecs in extent (Downes & Solomon 1998). Most of 
the luminosity of the galaxy is emitted by cool dust on larger scales. The 
equivalent radius of a blackbody sphere emitting 10 13 L at 50 K is 950 pc, 
several times larger than the observed sizes of nearby ULIRGs. 

The caustic curves predicted by the lens models that describe the image 
configurations observed in IRASF10214 (Broadhurst & Lehar 1995) and 
H1413 (Kneib et al. 1998) are about 0.8 and 0.2 arcsec in size in the plane of 
the source, larger than the scale of a 200-pc high-z source. Hence, the whole 
emitting region of the source should be subject to differential magnification; 
a difference in magnification by a factor of about 100 would be obtained 
between a 200-pc outer radius and a 2-pc inner radius. 

If the far-IR SED of a galaxy is modeled by a single-temperature dust 
spectrum, then differential magnification would be expected to increase the 
temperature for which the best fit was obtained, by a factor of about 10 
to 20 per cent for discrete data points at the wavelengths given in Fig. 1. 
This effect may account for at least part of the difference between the very 
high dust temperature of 107 K inferred for APM 08279 and the cooler dust 
temperatures generally inferred for other ULIRGs. 

The mid-IR SEDs of distant galaxies that are either unlensed or known 
not to suffer from differential magnification are very uncertain. This is likely 
to remain true for some years. In order to probe the spatial structure and 
temperature distribution of dust in the inner regions of AGN directly, a 
space-borne mid-IR interferometer with 500-m baselines at 20 /im will be 
required (Mather et al. 1998). In the meantime, it might be possible to use 
observations of variability to probe the same regions. The energy output 
from an AGN can vary strongly on short timescales. The response time 
of the hottest dust clouds, on sub-parsec scales, should be comparable to 
the light crossing time - about 1 yr. Over the lifetime of SIRTF, any such 
variations, amplified by differential magnification, should be detectable. 

4. Conclusions 

1. The slopes of the mid-IR SEDs of high-z lensed ULIRGs are likely 
to be flatter than those of unlensed ULIRGs because of differential 
magnification. While observations of these sources must be exploited in 
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order to investigate the nature of this class of galaxies, careful account 
must be taken of the potential uncertainties introduced by differential 
magnification, as well as by the effects of the orientation of a dust torus. 
Recent samples of submillimetre-selected galaxies should be immune 
from this problem. 

2. The SED of the longest known high-z ULIRG IRASF10214 is proba- 
bly affected by differential magnification in the same way as those of 
APM 08279 and H1413. However, it has a steeper mid-IR SED. This 
could be due to either a steep intrinsic SED or to a great optical depth, 
which obscures the most central regions even in the mid-IR waveband. 

3. The SEDs of a sample of lensed galaxies can be used to probe the 
conditions in the cores of distant dusty galaxies. A large sample of 
these objects will be compiled by the Planck Surveyor satellite (Blain 
1998). A direct test of the conditions in the central regions of ULIRGs 
will require a space-borne mid-IR interferometer (Mather et al. 1998). 

Acknowledgements 

I thank the staff of the Observatoire Midi-Pyrenees for their hospitality, 
and the MENRT for support while this work was completed. 

References 

Antonucci, R., 1993, ARA&A, 31, 473 

Barvainis, R., Antonucci, R., Hurt, T., Coleman, P., Reuter, H.-P., 1995, ApJ, 451, L9 

Blain, A. W., 1999, MNRAS, in press (astro-ph/9812415) 

Blain, A.W., 1998, MNRAS, 297, 511 

Blain, A.W., Longair, M.S., 1993, MNRAS, 265, L21 

Broadhurst, T., Lehar, J., 1995, ApJ, 459, L41 

Downes, D., Solomon, P.M., 1998, ApJ, 507, 615 

Eisenhardt, P.R., et al., 1996, ApJ, 461, 72 

Granato, G.L., Danese, L., Franceschini, A., 1996, ApJ, 460, Lll 
Green, S.M., Rowan- Robinson, M., 1996, MNRAS, 279, 884 
Hughes, D. et al., 1998, Nat, 394, 241 

Irwin, M. J., Ibata, R.A., Lewis, G.F., Totten, E. J., 1998, ApJ, 505, 529 
Ivison, R. J., et al., 1998, MNRAS, 298, 583 

Klaas, U., Haas, M., Heinrichsen, I., Schulz, B., 1997, A&A, 325, L21 
Kneib, J.-P., et al., 1998, A&A, 329, 827 

Lewis, G. F., Chapman, S. C, Ibata, R. A., Irwin, M. J., Totten, E. J., 1998, ApJ, 505, LI 
Mather, J. C. et al. 1998, astro-ph/9812454 
Peacock, J. A., 1982, MNRAS, 199, 987 

Rigopoulou, D., Lawrence, A., Rowan- Robinson, M., 1996, MNRAS, 278, 1049 
Rowan-Robinson, M. et al., 1991, Nat, 351, 719 
Sanders, D.B., Mirabel, I.F., 1996, ARA&A, 34, 749 

Soifer, B.T., Boehmer, L., Neugebauer, G., Sanders, D.B., 1989, AJ, 98, 766 
Soifer, B.T., et al., 1999, ApJ, in press (astro-ph/9810120) 
Xu, C. et al. 1998, ApJ, 508, 576 



